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I INTRODUCTION 



This report, is a theoretical and experimental analysis of a 
model (Wyngaard and LeMone, 1980) used to calculate the refractive 
index structure function parameter, > at the inter facial 
region at the top of an entraining, turbulent mixed layer. 

is related to the micrometeoro logical structure functions 
for temperature, humidity, Cq 2, and T-Q covariance, 

Cq»Q. The mixed layer inter facial region is important for EM 
propagation because is greatly increased by large T and Q 

fluctuations due to the entrainment of warm, dry air from the 
nonturbulent atmosphere above the mixed layer. 

Assuming that the rate of entrainment is in equilibrium with 
the free tropospheric virtual potential temperature (buoyancy) 
lapse rate, the model indicates that is proportional to 

(Ax) 2 9^* Zj_“2/3/ii0^ where X is T or Q, AX r.he jump 
at the interface, Zj_ the height of the ooundary layer and 

convective temperature mixed layer scaling parameter. 

The theoretical basis of this model is examined and four data sets 
from the NPS aircraft measurements program are used to test r.he 
mode 1 . 
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II THEORY 



A. Background 

The structure function parameters for temperature, Crp2^ 



and specific humidity, Cq^ , are to be evaluated in the 



inversion region by averaging between heights 2 = hg ana 
Z = h 2 (see Fig. i). The complete theory was developed by 
Wyngaard and LetMone (1980), hereafter referred to as WL, so only a 
summary of the derivation will be presented in this report. In a 
few instances WL’s work will be expanded to make certain 
assumptions and manipulations more explicit. 

The height h^ is defined as the top of the mixed layer 
where wd^ =0. At h 2 both fluxes and flux divergences are 
equal to zero. The average structure functions are 



where Ih = h 2 - Hq and the 0, 2 on the integral denotes Hq, h2 • 

The average structure functions are related to their 
respective dissipation rates by the Corrsin equation 




(la) 




(lb) 



< C-Tt — 1.0^—^ 1/'^<X^> 



(2a) 



= 1 , 0 ^— ^~ 1 ^ > 



'2b) 
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where e is the rate of dissipation of turbulent kinetic energy, X 
and Xq are the scalar dissipation rates (the factor 1.6 implies 
X is the rate of dissipation of temperature variance 9^). 

B. Evaluation of <X> 

For the moment, the development will be confined to the 
specific humidity (Q) . The dissipation rate is calculated from 
the scalar variance budget equation (Q denotes mean while q 
denotes fluctuating specific humidity; later in the paper q will 
denote mixing ratio, Q/p) . 

dv/dt + Wdv/dZ + d(wq)/dZ + 2 d(Q/p)/dZ = -XQ (3) 

where v = W is the mean vertical velocity (subsidence) and p 

is the density of air. Integrating this equation from h^ to 
h 2 > as in Eq . 1, yields 

<Xq> — -<Dq> - <Tq> - <Pq> (4) 

where D is the first two terms in Eq. 3, T (transport) the third 
and ? (gradient production) the fourth. Assuming "quasi-steady" 
conditions, WL show that <Dq> and <Tq> are negligible compared 
r.o <?,-^>; therefore 



<Xq> = -<?Q^ 



( 5 ) 



At this point the generalized inversion structure inodeL 



(Deardorff, 1979) is introduced 



Q = Qo + f(Z); 



ho<Z<h2 



(6a) 



dQ/dZ = Aq df/dZ 



(6b) 



where f(Z) describes the shape of the Q profile in the inversion 
region (assumed to be the same for C and T) with f(hQ) = 0 and 
f(h 2 ) = 1» Qq is the mixed layer value and ■iQ the jump in Q 
across the inversion. Substituting Eq. 6b into Eq. 5 and 
integrating by parts one obtains 



< P Q ^ Ah 



2 AQ 




d ( wq ) / dZ 



fdZ 



(7) 



The mean Q continuity equation 



-d(wq)/dZ = dQ/dt + W dQ/dZ 



is used in Sc . 7 to obtain 



r2 






-2 AC j iC/dt fdZ - 



2 Ac V zQ/ 



(9) 



The time derivative of Ea. 6a 



iC/dt = dO„/dr. * f dAC/ it 



( 10 ) 



J. 



and Eq. 6b can be subst.ir.ut.ed into £q. la. First the 
"quasi-steady" assumption is invoked, setting the following 
conditions 



dAQ/dt 


= 0 


(11a) 


dA0 y/ dt 


= 0 


(11b) 


dh^/dt 


= 0 


(11c) 


dAh/dt 


= 0 


(lid) 



However, since 

dho/dt = Wq + ‘^eo 



( 12 ) 



then Eq . 11c implies = -Wq . Assuming constant 

divergence 

W = Z/ho (13a) 



dW/dZ = W^/ho 



(13b) 



«2 = (1 + a) Wq (13c) 

where a = Ah/hc is the normalized thickness of the interfacial 
region. Employing these relations in Eq. 9 and doing the second 
integral by parts gives 

-<P^>lh = -21Q dc^/dt + ( AQ) -Wgc(l-'^- ^Q) ( 14 ) 



where the interracial functions Yo and 



are 



(15a) 



Yo = Ah 



Zq = Ah 



-^L 



f dZ 



f2 dZ 



(15b) 



The time derivative terra in Eq. 14 is eliminated by 
integrating the conservation equation (Eq. 10) from hQ to h2 



AhdQo/dt - Wqq AQ(1 + a - uYq) = wq, 



iO 



(16) 



which is substituted into Eq. 14 to obtain (WL Eq. 42) 
-<PQ>Ah = -2 AQYq + ( AQ) 2weoC-2YQ(l+a-aYQ) + (l + a-aZQ)] 



Later in their paper, WL use the equation 
-AQ Wqo (1 + a - uYq) = 



(18) 



which, in view of Eq. 14, obviously implies dQ^/dt = 0, Since 
WL have already required that dAQ/dt = 0, this solution appears to 
be quite restrictive. If Eq, 13 is used in Eq. 17 then 



w 



Ah 



(AQ)“ (1 






(19) 



Despite the simplicity of Eq , 19, WL prefer to keep the w-Gq 
term separate in their development. The primary reason for this 
is to simplify the analagcus development for since 
'^^vo ~ Therefore, WL now employ the "quasi-steady” 

entrainment formula 



= 0.8 W* S-l/{l+a) 



( 20 ) 



where 



W 



eo 



S = gF 92 h2/{w*2T.) 



with r 02 = dQ-^/dZ at Z = Y 12 snd W* is the convective 
scaling velocity (Zj_ = h^/ 0 . 8 ) 



W*3 = g ^ 0 ^^ 2^/T 



( 21 ) 



( 22 ) 



and "s" denotes the surface value. 

Rather than make an explicit substitution for ’^qq at this 
point, one could keep W^q as a variable, giving 

<Xq> = -2AQ Yq v^/(aho) + (AQ)^ Weo(l+a)SFQ/ho 

where 

Fq = C-2Yq (1+a-aYQ) + ( 1 +a-aZQ ) ]/ { a ( 1-^a )S ) 

Using the WL solutions obtained for Eq. 24 (and Yq = 1/2) 

Fq = (6R)"^ 

then 

<'<Q> = ~AQ wq^ YQ/'(ahQ) *■ ( AC ) ( 1 -t-ct ) S / ( oRh^ j 

Where R = g ) 

The final result, is obtained by substituting for •'^q in Eq. 
First the 3,^ continuity equations ar Uq and at h 2 are 



(25) 

(26) 



(23) 

(24) 



26 . 



combined with the hQ to h 2 integral form similar to Eq. 16 to 
produce the relation 



d A0^/dt + Wqq Ah(l+a)r92 “ ^ eo AOvU + a-aYg) = wB^o (28) 



Since WL assume dA0^/dt = wB^q = 0, 



A0^(l+a-aYQ) = rg 2 h(l+a) 



(29) 



Using Eq. 29 and Sq . 1 8 in Eq. 26, one obtains 



<Xq> 



(AQ)2 Wgo(l+a)s (l+6-l)/(hoR) 



(30) 



Note that the first terra in Eq. 26 (which was proportional to 
wqQ ) is six times as large as the second term (WL obtain 15/2 
for this ratio because they use two separate formulae for 
which differ by a factor of 4/5, i.e. 6*5/4 = 15/2). 

The development for temperature is parallel until the 6^ 
equation analagous to Eq. 26 is reached. Since wB^q = 0, tJhe 
final result is 



C. Structure Functions 

The final sr.ep in this process is to specify that < s> is one 



0,^ ) 2 ( 1 ) WgQ S/{6Rh^) 



(31a) 



<X<2> = 7( Aq) 2(i+a)w^Q 3/(5Rh^) 



(31b) 
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half the value typically found at under convective conditions 



<£>1/3 = (0. 2)1/3 V/* 2j_“l/3 

Assuming the "quasi-steady" entrainment rate, the structure 
functions become 

<Cq2> = 3.9(Aq) 2 9^*/(Zi2/3i.0^) 

<Ctv^> = O.5A0^ 9y*/z^2/3 

where 9^* = The virtual temperature structure 

function is related to the temperature structure function, 
0^2, by 



<C/p2> _ 2 <Cfp2>/ A0y 

where is the function given by WL. 

One point worth more discussion is the approximation 
Fq = (6R)”1 and the final forms of Eq. 31. Suppose the 
results of Eq. 19 were used and a different function defined 

< Xq > = (AG ) ( 1 -ra ) S C-Q/'hQ 

<Xg> = (AGv)2w2Q(l+a) 3 FQ/'h^ 

where Fq remains as per Eq. 19 but. 



(32) 

(33a) 

(33b) 

(34) 



(35a) 
(3 3b) 



(36) 



= (l-+a- aZQ)/(a (i-ta )s) 

Using Eq. 29 one can show 



f’o = - 2 Yq/R (37) 

Following the calculations WL have in their Appendix A, a , 
Fg and Gg are unique function of R/ S (providing the mixed 
layer gradient is zero) . Gg is considerably less variable than 
Fg . The following formula are reasonable approximations for 
0.KR/S<10 



a = 0.96 R/3 - 0.11 (R/S)i*5 
Gg = (1 + 0.064 '/r/S)R“1 
Fg = (1 0.28 s! R/S)R"l/6 



(33a) 



(33b) 

(33c) 



These formulae lead to slight modifications to the structure 
function equations 



<Cg2> = 3.3 .'A Q) (3j_2/ 

"'^Tv ^ —0.57 ( A0 Y ' V./ *0*7 /■ 3 ^ 



where Dm = 1 



0 . 22 




(39a) 

(39b) 



rs 



2 



can be written in 



The equations for Cq^ and C/t> 



various general forms (now dropping the bracket notation) 



Cx2/((AX)2dxEx) = 1.14 9^*/(A0 ^2^2/3 ) 
or, without substituting explicitly for and e 


(40) 


0.53(l+a) r 
Cx2/((AX)2dxEx) ° ASr--<M//' 

V £ 

where Dq = 1, Eq = 3, and 


(41) 


Et = Ti/A0^ 


(42) 


D. Discussion 





It is of interest to ponder the significance of the various 
"quasi-steady” assumptions (Eq. 11, 12, 13). Suppose we exhume 

the original conservation equation integrals from Deardorff's (1979) 
paper (his Eq. 13 and 21). Assuming only horizontal homogeneity 
and constant divergence, the general equations become 



Ih dQQ,dt + AhYQ dAQ/dt - iCLd - Yq) ^ -Q'^eo- = '■dc 


(43a) 


Ah de,,^/dt ^ AhYQ dA0,,,/dt - Ae.,L(l - Yq) ^ '^Q •'''eo^ = ^ 


(43b) 



(43b) 



Thus Eq 16 and the 9v analogue can be reproduced by requiring 



Wg2 = (1 + ct) '^eo 

dAo/dt = dA0^/dt = 0 (44b) 

It is not necessary to require = -w, dhQ/dt = dh 2 /dt =dAh/dt = 0. 
This explains why WL found excellent agreement with Aschurch data 
where W = 0 and - 10 cm/s. 

Similarly, the general forms for the dissipation integrals 

are 

-<PQ>lh = -2 Ah AQYQdQ^/dt - 2 AhZQACdAQ/dt 

+ (AQ)2[{i - Zq) Wg2 + ZqWso^ 

-<P^>Ah = 2AhA0^ (Yq2 - Zq) dA0,^ /dt 

+ ( A0v)2C-2YqL(1 - Yq) Wg2 + Yq (45b) 

+ (1 - 2q) Wg2 Zq 

which reduce to Eq 14 if the conditions of Eq 44 are met. 

Since entrainment and surface flux tend to counteract each 
other in the Q case, it seems quite reasonable to assume that the 
dAQ/dt and dQ^/dt terms are negligible in Eq 45a 

-<?Q>Ah = (AQ)2 -(1 - Zq) Ws2 ^ 2q (46) 

Inst-aad of .naking r.he assumption Eq 44a, suppose ve simply assume 



awe vsMQv 



(47) 



We2 = Weo = 

which is the standard cloud-free result from Lilly (1968) where 
typically a =0.2. Then one can easily show that 

<Cq2> = 3.3 (aq) 2 9.^*/(Zi2/3^0^) 

which is identical to the WL result as expressed in Eq. 39ai In 
other words / the combination of "quasi-steady" assumptions 
We 2 = (1 + and Wq 2 =^vs/^”0^o^ 

equivalent to the assumptions of Eq. 47 even though they may imply 
vastly different entrainment rates. 

If one uses the assumptions of Eq. 47 and parallels the WL 
development, then the equivalent to Eq. 18 is 



'^eo (49) 



and the equilibrium condition form the 9 ^ equation is 



w0 yo “ ^"^eo 0 ” (50) 



which, assuming w0 = 0, gives 



The results for 9,, are also incerescing oecause it is not 
clear that the dAG,^/dt term should he negligible compared tc the 



other terms in Eq 46 b. Suppose we Let 

-<?0>Ah = A + B (52) 

Then the dA0^/dt term is small if A/B is small (returning to the 
“quasi-steady" format) 



A/B 



ho(YQ2 - Zq) dA0^/dt 6(R/S) 
A0V (1 + a) Wgo 



Since Yq2 - Zq 



-0.1, we can write 



A/B 



-0.6 hQ(R/S) dA0^/dt 

(1 + a) Wqq A0^ 



(53) 



(54) 



The magnitude of A/B can be examined by using the general 
relation 



dA0^/dt = -d0Q/'dt -r 



(55) 



and writing a simple entrainment formula (e.g. "quasi-steady" 



•'< o 2 ~ V s ' ^ • A ‘^o ^ 



V 30 j 



The inr.ecral of -he conservation ecuation from 3 = 0 to 






ves 



therefore 



h(3 d0^Q/dt — w0^g 






(57) 



d/10y. /dt = - WgQA0v/hQ (58) 

using Eq. 54 we find 

0.6 

A/B = — R/S (59) 

( 1+a ) 

A good example is the Aschurch data quoted by WL where Eq. 57 
was shown to be applicable. Since R/S =0.3 for that data, A/B = 
0.15 and dA0y/dt is negligible. 

Certainly the conditions set by VJL are consistent with 
neglecting dAO^/dt. It is not clear how to identify conditions 
where this assumption is invalid. Eq . 54 cannot provide much 
guidance because it is based on solutions to Eq. 23 with 
dA0^/dt =0. It is interesting that in the conditions where the 
WL equations for "quasi-steady" entrainment are expected to 
breakdown (A0,/ large, R/S >1) then the Lilly type relations give 
the same results for Cq 2 . xf dAOy/dt terra becomes 

important, then one anticipates the WL formulation will 



underesti.mate Co.2. 



Ill ATMOSPHERIC DATA 



A. Measurement Techniques 

The measurements were made using a single engine Bellanca 
Viking aircraft operated by Airborne Research Associates of 
v^eston, MA. The instrumentation and data processing have been 
previously described in detail (Fairall et. al^ 1980; Schacher 
et. al., 1980) so only a brief summary is given here, 

i) Mean temperature, T: platinum resistance sensor with 

standard aircraft mount. 

ii) Mean humidity, Q: cooled mirror dew cell. 

iii) Mean windspeed, U: estimated at the surface from the 

sea state and DMV navigational aid. The present LORAIJ system 
was not available. 

iv) Sea surface temperature, TgZ Barnes PRT-5 IR 
radiometer . 

v) Cm2; microthermal sensors (4.5 um dia. tungsten) 
in the paired configuration. 

vi) Cq 2; Lyman-alpha fast humidiometers using the 
inertial subrange filter method. Absolute calibration based 
on comparison with a microwave ref ractometer . 

vii) £: hot wire (4.5 um dia. tungsten) constant 

temperature anemometer. The inertial subrange filter method 
was used. 

3. Surface Fluxes and Turbulence Scaling Parameters 

Surface fluxes were evaluated from aircraft measurements 
using two methods: a) bulk aerodynamic and b) dissipation 

(inertial subrange). The fluxes are defined in terms of the 



lb 



following scaling parameters; 



momentum: 



p uWg = -pu*2 



(60a) 



sensible heat; pCp w0g = -pCp u*T* 



(50b) 



moisture; 



oqWg = -pu*q* 



(50c) 



The momentum flux is also referred to as the Reynolds stress, T. 
Note; the bul)c method was not used overland. 

1. Bul>: aerodynamic method. 

The exact details were described in a recent paper (Davidson 
et 1981). Using Eq. 4a from that paper, one can relate the 

values of some meteorological variable (temperature, moisture or 
wind speed) at the sea surfce, Xg, and at some height Z in the 
surface layer, X^, to the scaling parameter, X*; 



where 3 q and Zq-j are roughness lengths, L is the Monin-Cbukhov 
length, 3 and k are constants, and and yr], are empirical 
functions . 



“ u 2 C in ( Z / Z Q ) 



(2/L)]-l 



(61a) 



T* = (T, - Tg) SkCm (Z/Zq.t) - ’^'t 



(61b) 




(51c) 



2 . 



Dissipation method. 



The dissipation method relies on semi-erapiricaL relationships 
of inertial subrange turbulence to surface-layer scaling 
parameters (Fairall et ai. , 1980). The equations are 



u* = C( £ k Z)/4i (Z/L)]1/3 
T* = [z2/3 C^^/f { Z/L) ]^/2 
Q* = Cz2/3 Cq2/(A f(2/L)]l/2 



(62a) 

(62b) 

(62c) 



where £ is the dissipation rate, (p and f are empirical functions, 
and A. is a constant. Because the structure- function parameters 
Crp2 and Cq 2 are related to the square of the scaling 
parameter, a sign ambiguity exists. This can usually be 
eliminated by examining the low-level height dependence of s, 

Cq 2 and Crp2 because the functions ^ and f have 

character istic profiles for stable and unstable conditions. 

Both methods yield accuracies on the order of 10% for u* , 
jr0.02®C for and ;^0.02 g/m^ for Q* (note: q* = Q*/o). 

C. Bata Sets 

The data given in this report were cctained in four field 
programs: 

i) Panama City (PC), 1978 (more detail available in 

Fairall, 1979) over the Gulf of Mexico in Florida. 

ii) Wliite Sands (WS), 1979. Two profiles over the desert 

under highly convective daytime conditions. 
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iii) MAGAT (MG), 1980 (more detail available in Fairall, 
1980) in the Monterey Bay area. 

iv) Bahamas (BH), 1980. A series of profiles taken near 
Andros Island in the late fall. 

The complete data sets were examined to remove profiles tliat 
encountered boundary-layer clouds. A total of 23 profiles were 
selected. Graphs of the mean and turbulence profiles for each 
case are given in Appendix A. A summary of the basic scaling 
parameters is given in Table 1. 
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TABLE 1. 



Meteorological data and surface scaling parameters 
from the cloud free NFS data sets. 
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IV. RESULTS 



A summary of the secondary scaling parameters used for the 
NFS data set is given in Table 2. Also shown in Table 2 is a 
comparison of the measured and model assumed values for z at the 
inversion. With very few exceptions, the model assumption (Eq. 
32) is very good. The entrainment velocities calculated from the 
''quasi-steady'* assumption used by WL (Eq. 20) and the more 
conventional parameterization of Lilly (1968). 

W^o/W* = 0.2 S.^*/A0^ (63) 



are also calculated. 

In Table 3 are the measured values of C<pA and Cq 2 at 
the inversion plus their normalized forms 

Ix = Zi2/3 Cx2/((Ax)2 Dx ?x) (64) 

taken from Eq. 40. According to WL (Eq. 26), the theoretical 
value is 

I^ = 1.14 3.^*/A0.^ (65) 

which is the same for T and Q. 

.A direct comparison of .measured and calculated values of 
C.-p2 and Cq 2 given i.m ?ig. 2. The medal pradiems ohe 
measure.ments within a faemor of three. The uncertainly is 
slightly greacer than the factor of two suggested by WL but 



Table 2. 



Surface scaling C^vs cori^ective scaling 

(W*, 9^* and and inversion scaling (R, S 

and Wqq) parameters. Two formulae are used to 
estimate "steady" is Eq 20 and "Lilly" is Eq 

64. 
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Table 3. 



Measured values of the interfacial structure functions 
(C<p2 and Cq^) and their resultant values for 

^x^where X=^ • 

or Q Tnese are compared v/ith theoretical values, 1^,, 
using the "steady” and "Lilly" entrainment values. 
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includes various measurement, errors and uncertainties. 



I'lote that 



the Cq 2 data has a greater range of values than . This 
is consistent with the WL model. If we examine the function 

H = Zi2/3 Cx2/(Dx 9 V*) (66) 

then 

Ht = 10 V (67a) 

Hq = (AQ) 2/A0 ^ (67b) 

A graph of Hip and Hq is shown in Fig. 3 for a typical range of 
A0V and AQ from the NFS data set. Note that H-p varies roughly 
from 2 to 9 while Hq varies from 4 to 72. 

The entrainment parameterization was tested (Fig. 4) by 
plotting measured values of Ix(H;q. 65) against the model value 
( Eq . 66) which is based on the entrainment formula given by WL 
(Eq. 26). This plot gives a much higher correlation than a 
si.milar graph (not shown) using the more traditional formula due 
to Lilly (1968), Eq. 62, which gives 



(Lilly) = 0. 



(1 



■ 4 



(A0’/^^ 



(63) 



This is not really significant because, when used in proper com- 
bination with Eq. 48, the Lilly formulation also leads to Eq. 66. 

In order to look for systematic errors, the ratios ( Rrp and 
3.Q ) of .measured to .model values of Cp2 ^nd Cq 2 vere 
calculated and plotted against .10 y (Fig. 5). si.mp le 
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log-average yields R.p = 1.15 and Rq = 
suggest.s that, the model underestimates C<j.2 (large R<j< ) when 
A0y is small while it overestimates when AOy is large (the 
Cq 2 data is too sparse to clear up this question) . This could 
be due to an error in the estimation of A0^ and AQ (admittedly 
rather subjective) . An e.xamination of Fig. 3 suggests that a 
reasonable adjustment of A0y (several tenths K) will not move 
the data points substantially closer to the Rip = 1 midline. 
Another possibility is that Eq. 20 tends to overestimate 
when A0V is large while underestimating for small A0y» 

Given the considerable scatter in the results, the 
uncertainties in the esti.mation of A0v snd AQ from measured 
profiles and the insensitivity of C<t. 2 ^^id AQ it is 

suggested that a simplified formula for can be used for 

application to radiosonde quality data. If one assumes (based in 
Fig. 3) that Hrp ~ 5, then 

Cp.2 :: 5 3^* Zi~^/^ (59) 



Based on 
leasr. as 



nhe NPS data set this approximation appears 
accurate as the niore complicated formula {? 



to 



9 • 



6 ) 



a t 



AGvlKl 



Rj (Meas./ Theory) 



o 




? I jure 6. The ineasure'/i value of iivi-ied oy rne mcoel 

vaLue usinc ^he siioplified ercpression (Eq. 70). 




V CONCLUSIONS 



The Wyngaard-LeMone inversion Layer scaling has been examined 
theoretically and tested against a data set obtained by NPS 
investigators in cooperation with Airborne Research Associates. 

The theoretical examination indicated the following: 

i) The WL theory is more general than is implied by the 
strict assumptions of the "quasi-steady" theory. 

ii) The WL development can oe simpLifiec slightly, leading 
to modest adjustments of the normalization constants. 

iii) The steady state assumption that dAQ/dt is negligible is 
reasonable under most conditions. The assumption that dA0^/dt 

is negligible may not be justified when R/S > 1. 

The examination of the atmospheric data indicated the 
fo llowing: 

i) The assumption that e at the inversion is proportional 
to a fixed fraction of the surface buoyancy flux was quite 
reasonable . 

ii) The WL model predicted the measured value of C,p2 ^^d 
Cq 2 to within a factor of three. 

iii) Some evidence, though statistically weak, was found to 
suggest the model overestimates the structure functions for Large 

( > 3K; ) while it underestimates for small ■. < 2I<) . In 

the other hand, this could be a manifestation of the Stei.t effect 
for comparison of data sets subject to error vhere small values 
are usually overestimated and large quantities are usually 



underestimated . 



Based on these results, it is obvious that a raajor weakness 
of the raodel is its reliance on an entrainment formulation that is 
too restrictive. The two extremes of the buoyancy jump (A0v^ 
may involve different entrainment regimes (e.g. encroachment, 
convective instability or the Lilly formulation). It would also 
be useful to include the effect of inversion windshear on Wg and 
on the structure functions . Another area of investigation might 
be stable surface layers. These may be very important for surface 
optical propagation because C<p2 values are often sizeaole and 
2j_ is usually small (on the order of 100m) . 
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APPENDIX A 



This appendix contains graphs of mean (©y/ q) and 
turbulence (Crp2^ £) profiles for each of 23 data 

sets. The site designations are defined in Section III-C. The 
abstraction of 
1, 2 , 3 in the 



this data to obtain the relevant parameters (Tables 
main next) is described in Section III. 
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